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ABSTRACT 

We investigate the density profiles and evolution of weakly self-interacting cold dark 
matter haloes using a numerical code based on the collisional Boltzmann equation. 
This approach is alternative to N-body techniques in following the dynamical evolution 
of haloes in the cosmological context and taking into account particle self-interaction. 
The physical case with a cross section inversely proportional to the relative velocity 
of the colliding particles is modelled with an unprecedented resolution, spanning five 
orders of magnitude on the radius for each halo. The modelled haloes cover a mass 
range from dwarf galaxies to galaxy clusters. We find that for cuioo ~ 10~ 24 cm 2 /GeV, 
where a is the cross section per unit mass and Vioo is the collision velocity in units 
of 100 km/s, soft cores in good agreement with observations on galactic as well as on 
galaxy cluster scales are obtained. Remarkably, the observed nearly invariance of the 
halo central density with mass is reproduced. 
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1 INTRODUCTION 

The current Cosmology based on the growth of small fluc- 
tuations through gravitational instability in a universe with 
cold dark matter (CDM) defines a physically well-motivated 
theory, capable to explain most of the properties of the 
large-scale structures of the universe. In this scenario, where 
a large fraction of the matter is supposed non-dissipative, 
cold and collisionless, the dark haloes form hierarchically 
via gravitational collapse of primordial density fluctuations. 
The density profiles of these dark haloes have been studied 
systematically using N-body simulations (Navarro, Frenk & 
1996, 1997; hereafter NFW). It was found that a simple an- 
alytical density profile with one free parameter is able to 
describe the typical structure of the virialized CDM haloes. 
On its own, this free parameter, which can be the concen- 
tration, depends only on the halo mass. 

The NFW density profile goes as p oc r _1 at small 
radii. According to recent high-resolution simulations, the 
inner density profile of the haloes is even steeper than r -1 
(p oc r" 1 ' 5 ) (Moore et al. 1999; Jing & Suto 2000). This 
prediction of the CDM scenario is in potential conflict with 
observations: the halo inner density profiles inferred from the 
HI rotation curves of dwarf and low surface brightness (LSB) 



galaxies and from a mass map constructed by strong lensing 
techniques for the galaxy cluster CL0024+1654 seem to be 
much shallower than r _1 (Moore 1994; Flores & Primack 
1995; Burkert 1995; de Blok & McGaugh 1997; Tyson et 
al. 1998) . Recently, the rotation curves of some of the dwarf 
and LSB galaxies were re-obtained using high-resolution Ha 
observations, resulting the inner density profiles more con- 
centrated than those obtained from the HI observations, but 
on average still shallower than r _1 and with a large scatter 
(Swaters et al. 2000; van den Bosch & Swaters 2000; Dal- 
canton & Bernstein 2000). Besides, as Firmani et al. 2000b 
have shown, the halo inner density profiles of most of the ob- 
served LSB galaxies could be shallower than reported if the 
baryon contraction is taken into account (see also Firmani 
& Avila- Reese 2000). 

An interesting property to be confirmed with more ob- 
servational studies, in particular for galaxy clusters, is that 
the central density of dark matter haloes is roughly inde- 
pendent of the halo mass with a value around 0.02 Mq pc -3 
(Firmani et al. 2000a,b). Not only the predictions of the 
CDM scenario, but also of alternative scenarios as the warm 
dark matter one are unable to predict such a scale invari- 
ance of the central density of dark haloes (Avila- Reese et al. 
1998; Coh'n et al. 2000). 
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In a recent burst of activity, several alternatives to the 
well motivated CDM scenario were proposed in order to 
overcome the difficulties mentioned above and others as the 
excess of satellites and a little angular momentum in the 
galaxy disc predicted by numerical simulations. These alter- 
native theories try to suppress the excess of central mass of 
the galactic haloes and may be summarized into two classes: 
1) a manipulation of the power spectrum with a suppression 
of the small scales as an intrinsic property of a dark mat- 
ter: Avila-Reese et al. 1998; Moore et al. 1999; Firmani et 
al. 2000a,b; Peebles 2000; Colin et al. 2000. These works 
make different assumptions concerning thermal energy at 
the beginning of the mass aggregation history (MAH). 2) 
Assuming for dark matter particles a different nature from 
collisionless: self-interacting (Spergel & Steinhardt 2000); re- 
pulsive (Goodman 2000); decaying (Cen 2000); with a lim- 
iting phase-space density (Hogan & Dalcanton 2000) and 
annihilating (Kaplinghat et al. 2000). Intriguingly, only the 
assumption of a self-interacting or annihilating cold dark 
matter seems to be able to reproduce the observed central 
density scale invariance ranging from galaxy to galaxy clus- 
ter size. 

In a previous work, adopting the self-interacting CDM 
scenario suggested by Spergel & Steinhardt, we have pro- 
posed a thermodynamical model able to describe the 
gravothermal expansion of the core in the haloes (Firmani 
et al. 2000a); with this model, the scale invariance of the 
halo central densities was predicted. The NFW density pro- 
file shows a velocity dispersion raising outwards with radius 
(Cole & Lacey 1996; Fukushige & Makino 1997). Therefore, 
if particles are self-interacting, the inward heat transfer in- 
duces a thermalization process able to avoid the formation 
of a cuspy core. The influence of weak self-interaction on the 
structure of the CDM haloes has been recently explored us- 
ing cosmological numerical simulations (Yoshida et al. 2000; 
Dave et al. 2000). These works, where constant cross sec- 
tions were used, have led to interesting results, that we will 
discuss in comparison with ours, later. 

A cross section inversely proportional to the particle 
velocity dispersion seems to be better motivated by obser- 
vations (Firmani et al. 2000a,b; Wyithe et al. 2000). In this 
paper we focus on detailed dynamical predictions of the den- 
sity profiles of self-interacting spherical symmetric dark mat- 
ter haloes ranging from galaxy to galaxy cluster mass. The 
cosmological model we assume is flat with vacuum density 
Qa = 0.7 and expansion rate Ho = 65 km s _1 Mpc -1 . 



2 THE COSMOLOGICAL BOLTZMANN CODE 

Our dynamical approach aimed to investigate the cosmo- 
logical halo density profiles is based on a solution of the 
collisional Boltzmann equation for dark matter. The viri- 
alized dark halo is described by the distribution function 
N(E, J), where N(E, J) ASA J gives the number of parti- 
cles within the intervals AE and AJ centred at the total 
energy E and the angular momentum J. The particle popula- 
tion is assumed to be uniform with respect to the directions 
of J. The mass aggregation supplies new particles, each one 
with a total energy and angular momentum. The computa- 
tion of each particle orbit is based on a gravitational field 
derived from the mass distribution in a self-consistent way. 



The change in time of the gravitational field due to the mass 
aggregation modifies the particle orbit in agreement with its 
adiabatic invariants. The particle mass is distributed along 
its orbit according to the fraction of the period spent in any 
interval. For collisionless particles this approach is similar 
to the secondary infall (Avila-Reese et al. 1998) and allows 
to calculate a virialized halo density profile for mass and 
angular momentum aggregation histories in a specific CDM 
universe. For a given universe, the mass aggregation his- 
tory (MAH) is obtained making use of the extended Press- 
Schechter formalism based on the conditional probabilities 
for a Gaussian random field. For a given mass, we generate 
a set of MAHs through Monte Carlo simulations and we cal- 
culate the average MAH. If the perihelion to aphelion ratio 
of the orbits (angular momentum) is opportunely selected, 
the NFW profile is obtained on galaxy and galaxy cluster 
scales. 

The effects of the collisions between particles is de- 
scribed by the collisional term of the Boltzmann equation 
and induce a further evolution on the distribution function. 
Our scheme will provide an accurate solution only when the 
collision time is sufficiently larger than the orbital period 
of the particles. The collision cross section multiplied the 
relative velocity of the colliding particles is assumed to be 
constant in agreement with the previous discussion in Fir- 
mani et al. 2000b. This hypothesis makes particularly ef- 
ficient the Monte Carlo method implemented to describe 
the collisions of the particles of each discrete state (E,J). By 
each trial we move a large amount of particles from one state 
to another according to the probability of the process. The 
results show an enormous space resolution obtained with 
a minimum computational expense, with each run taking 
approximately a couple of hours on a 450MHz pentium pro- 
cessor. 

The efficiency of our method is counterbalanced by the 
practical difficulty to explore solutions with collision times 
lesser than the particle orbital periods (high cross section), 
as well as halo structures which deviate from the spherical 
symmetry. For these cases complementary studies with N- 
body techniques have been planned. 



3 SELF-INTERACTING DARK MATTER 
HALO EVOLUTION 

We have calculated the evolution of haloes in a ACDM cos- 
mology assuming a self-interaction cross section per unit 
mass (a) given by: a vioo = 10~ 24 cm 2 /GeV, where vioo 
is the relative velocity of the colliding particles in units of 
100 km s~ . Fig. 1 shows the present-day density profiles of 
haloes with masses of 6 10 11 and 6 10 15 Mq (solid lines). 
For comparison, the corresponding NFW profiles are also 
shown (dotted lines) . Remarkably, the self-interacting CDM 
haloes appear to be thermalized in the centre resulting in a 
non-singular isothermal sphere, while in the outer parts, the 
density profile is well described by the NFW shape which is 
representative of the hierarchical merger history. The non- 
singular isothermal profile obtained with our code confirms 
the halo profile (dashed line) we derived in a previous work 
in order to describe observations (Firmani et al. 2000b). In 
that case, the hydrostatic equilibrium was assumed in or- 
der to pass from a NFW profile for the outer regions to an 



© 0000 RAS, MNRAS 000, 000-000 



Modelling Self- Interacting CDM Haloes with a Cosmological Boltzmann Code 3 




Figure 1. The behaviour of the mass density profile for haloes 
with total mass: M = 6 10 11 and 6 10 15 Mq (solid lines) with 
the self-interacting cross section: a fioo = 10 — 24 cm 2 GeV -1 . 
Dotted lines show the NFW profile for both haloes. Long-dashed 
lines draw the non-singular isothermal profile obtained in Firmani 
ct al. 2000b (sec the text) 

isothermal sphere in the centre. For that "naive" model the 
core-radius was obtained once assigned the central density. 

With the dynamical model based on the solution of the 
transfer equation we presented here, we can now estimate 
the typical number of collisions in the halo core over a Hub- 
ble time (i//): N co u = pavtn- For a 6 10 M halo and 
the assumed cross section, we find N co n « 4. 

Fig. 2 shows the evolution of the density profile of an 
halo of 6 lO n M0. Since early epochs (z ~ 40) a soft core 
of few pc size is formed. At z — 20 the central density is 
about 4 Mq pc~ 3 and the core size is already of ~ 40 pc. 
The central density is high enough to allow collisions on time 
scales of only a few hundred million years. Due to the mass 
aggregation, the inner radial velocity dispersion gradient is 
preserved. Then, due to heat transfers inwards, the soft core 
continues growing while the central density decreases. At 
any time, the core density is determined by a self-regulating 
mechanism which makes the collision efficiency able to trans- 
fer enough energy to expand the core. The central density 
decreases and the core size rises, both reaching values com- 
parable with the present observations, i.e., 0.03 Mq pc~ 3 
and 6 kpc, respectively. 

Of great interest is the invariance of the central den- 
sity, shown in Fig.l, in passing from an halo of galactic size 
to an halo with mass four order of magnitude larger, repre- 
sentative of a cluster of galaxies. In Fig. 3a and Fig.3b we 
compare the central densities and core radii predicted by 
the models with those inferred from observations of dwarfs 
(filled squares), LSB galaxies (empty squares) and the clus- 
ter CL0024+1654 (filled circle) (see Firmani et al. 2000b for 
details). Big empty squares are for those LSB galaxies with 



Figure 2. The evolution in rcdshift of the density profile is shown 
for an halo with M = 6 10 11 Mq. 

Ha rotation curves (Swaters et al. 2000). The models were 
calculated for two cross sections: crvioo = 10~ 24 and 2 10~ 24 
cm 2 GeV -1 (dashed and dotted lines, respectively). 

The agreement of the self-interacting CDM model with 
observations is rather remarkable. In a range of scales from 
dwarfs to galaxy clusters, the haloes have (1) a nearly con- 
stant central density, and (2) a core radius roughly pro- 
portional to the maximum circular velocity. We stress that 
the only way to obtain halo central densities independent of 
their masses is assuming a oc 1/v. 

It is interesting to note that constraining the density 
for the cluster to the value of p c = 0.02 Mq pc -3 as inferred 
by the observation of CL0024+1654, we derive an estimate 
of the cross section: ct«ioo ~ 10~ 24 cm 2 /GeV. 



4 COMPARISON WITH PREVIOUS WORKS 

In Fig. 4 we show a grid of models for two different masses: 
6 1O U M and 6 1O 15 M and wide range of cross sections. 
In order to scale the models with the mass for a given value 
of avioo, Fig. 3 shows that p c is roughly scale invariant while 
R c is proportional to V m . 

With the aim to compare our models, derived with the 
assumption a oc v^ 1 to models with a — const, we will as- 
sume the latter velocity dispersion as representative of the 
collisional velocity. From Fig. 2 of Yoshida and co-workers 
(2000), extrapolating the profile towards the center for the 
model SlW-a, with the same pattern of SlW-b and SlW-c, 
we read a central density of p c = 2 - 3 lO 16 M /h 2 /Mpc 3 
and a core radius R c = 35 — 50 kpc/h for a halo of 
7.4 10 14 M /h. We have simulated a cluster of M = 
6 10 15 Mq. Rescaling the mass of the Yoshida et al. (2000) 
halo (7.4 10 14 M Q /h and v = 1500 km/s) to the cluster mass 
used here (6 1O 15 M and v = 2400 km/s), it is easy to ob- 
tain p c = 0.013 - O.OO8M /pc 3 and R c = 85 - 120 kpc. In 



© 0000 RAS, MNRAS 000, 000-000 



4 Firmani, D'Onghia & Chincarini 
o 



o 

Q_ 



O 



-1 - 



o 

Q_ 
J*. 



□ □ 



1.5 



2.0 



2.5 3.0 
Log V m [km/s] 



3.5 



0.0 



" -1.0 



-2.0 - 



-3.0 



2.0 



1.0 



0.0 



-1.0 



6 10 11 M,-. 



6 10 15 M,-. 



6 10 15 IVL 



6 10 11 M,. 



-26 



-25 -24 -23 

Logov 100 [cm 2 Gev 1 ] 



Figure 3. The plot shows a comparison between central densi- 
ties and core-radii predicted by the model and the same quantities 
inferred by available observations of dwarfs (filled squares), LSB 
galaxies (empty squares) and the cluster CL0024+1654 (filled cir- 
cle) (see Firmani et al. 2000b for details). Big empty squares rep- 
resent LSBs Ha rotation curves (Swaters et al. 2000). Models are 
shown for a v 100 = 10~ 24 (dotted lines) and 2 10~ 24 cm 2 GeV" 1 
(dashed lines). 



our model, for the cross section of SlW-a, a — 0.1 cm 2 /g, 
we find p c = 0.01M Q /pc 3 and R c — 100 kpc. Let us now 
analyze the case SlW-b of Yoshida et al. (2000). In this 
model, which resolution is higher than the SlW-a case, the 
central density is p c — 4 10 15 M Q h 2 /Mpc 3 and the core ra- 
dius R c = 150 kpc/h. Making the same mass rescaling we 
obtain p c = 1.7 10~ 3 M o /pc 3 and R c = 370 kpc. This case 
is for a cross section a — 1 cm 2 /g. In Fig. 4 we see that 
our model for the same a predicts p c = 2 10 _3 Mq/pc 3 and 
R c = 355 kpc. Taking into account the different approaches 
the agreement is really outstanding. 

Regarding Dave et al (2000) work, the only case we 
have simulated is the first top panel from left to right of 
Fig.l of the same paper. For a — 10~ 24 cm 2 GeV~ x and for 
an halo of M = 6.17 IO^Mq they find p c = 0.03M Q /pc 3 in 
agreement with our value of central density of 0.03M Q /pc 3 . 



Figure 4. The top panel shows the central densities of virializcd 
haloes of: 6 10 11 M Q (dotted-line) and 6 10 15 M Q (dashed-line) as 
a function of the strenght of the self-interacting cross section. In 
the bottom panel core radii are plotted for the same haloes. 



5 SUMMARY AND CONCLUSIONS 

We have calculated the evolution of the density profiles 
of self-interacting CDM haloes using a numerical approach 
based on the collisional Boltzmann equation. This study 
presents two original points: (i) we were able to explore 
the halo density profiles covering about five orders of mag- 
nitude in the radius and for haloes from galaxy to galaxy 
cluster scales; (if) we used a self-interacting cross section in- 
versely proportional to the collision velocity of the particles. 
Analysing the results obtained with our code we conclude: 

(1) The interior of dark haloes may be deeply affected by 
self-interaction. A modest cross section value (avi m w 10~ 24 
cm 2 GcV -1 ) is enough to produce soft cores. 

(2) Collisions between dark particles induce a thermaliza- 
tion process with the consequent appearance of a soft core. 
Thermal equilibrium in the core is established with a few 
collisions in a Hubble time (« 4). 

(3) Soft cores appear already at the very early dynamical 
evolution (z > 40); with time the central density decreases 
and the core size increases. 

(4) The observed scale invariance of the central density 
from galactic to galaxy cluster scales is reproduced only if 
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the cross section is inversely proportional to the collision 
velocity. 

(5) The core radius predicted by the model increases pro- 
portional to the halo maximum rotation velocity in agree- 
ment with the observed trend. 

(6) The observed central density of CL0024+1654 is con- 
sistent with assuming a self-interacting dark matter cross 
section given by: ctdioo ~ 10 -24 cm 2 GeV -1 . 

The code is limited in analysing dark halo properties which 
depart from spherical symmetry. Further studies concerning 
dark halo triaxiality and the evolution of substructures have 
been planned using N-body techniques. 
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